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Abstract. Using the delta correction to the standard free energy [4] in the 
elastic setting with a quadratic foundation term and some parameters, we in- 
troduce a one dimension only model for strong segregation in diblock copoly- 
mers, whose sharp interface periodic microstructure is consistent with experi- 
ment in low temperatures. The Green's function pattern forming nonlocality 
is the same as in the Ohta-Kawasaki model. Thus we complete the statement 
in [31, p. 349]*: "The detailed analysis of this model will be given elsewhere. 
Our preliminary results indicate that the new model exhibits periodic mini- 
mizers with sharp interfaces." We stress that the result is unexpected, as the 
functional is not well posed, moreover the instabilities in L 2 typically occur 
only along continuous nondiffcrcntiablc "hairs" . 

We also improve the derivation done by van der Waals and use it and the 
above to show the existence of a phase transition with Maxwell's equal area 
rule. However, this model does not predict the universal critical surface tension 
exponent, conjectured to be 11/9. Actually, the range (1.2,1.36) has been 
reported in experiments [21, p. 360]. By simply taking a constant kernel, this 
exponent is 2. This is the experimentally (±0.1) verified tricritical exponent, 
found e.g., at the consolute 0.9 K point in mixtures of 3 He and 4 He. Thus 
there is a third unseen phase at the phase transition point. 



Key to some recent developments in modern microelectronics has been the ability 
to create high-quality atomically abrupt interfaces between different semiconduc- 
tors, which can produce new quantum states and uncover unexpected phenomena 
(see http://en.wikipedia.org/wiki/Nanotechnology). Several experiments contra- 
dicting common knowledge about bulk materials having only diffuse interfaces have 
been reported, see the References, especially [7], where a new experimental tech- 
nique is explained. 

It is remarkable that chemically diverse diblock copolymers and microemulsions 
admit very similar characteristic structures, e.g., the bicontinuous ordered double 
diamond [27]. A theory formally predicting the phase transition between weak and 
strong segregation was derived in [4] . Here we show how to use elasticity and choose 
parameters in the Vaserstein pseudoassociation potential to get local L 2 minimizers 
resembling the square wave [38] . As a by-product we get a rigorous proof of the 
gas to liquid phase transition. 

To make this note self-consistent, we review the history of this modeling. 

The ideal gas law 

(CD P-f. 



* A preliminary version of this result was announced at the "Defects and their dynamics" 
workshop (09.08-16.08.2003), Banff International Research Station, Alberta, Canada. 

1 



2 



ADAM CHMAJ 



where P is the internal pressure, V the volume per mole, T the absolute tempera- 
ture, R the gas constant, does not predict a phase transition between gas and liquid, 
defined according to experiment as two densities coexisting at the same pressure (see 
e.g., the isotherms for carbon dioxide in [36]). 

Van der Waals formulated the modified equation of state as 

(0 2) P - — - — 

[U - Z) V-b V 2 ' 

where a is the attraction parameter which arises from polarization of molecules into 
dipolcs and b is the volume enclosed within a particle (the repulsion parameter). 
Below the critical temperature isotherms for (0.2) have a wiggle [36], which is 
unphysical. Maxwell's construction in which the fluid is taken around a reversible 
cycle of states from a logical point of view is worthless [36] , as states on the wiggle 
have no meaning. However, such soft reasoning can give other interesting results 
[22, 36]. 

For uniform systems, the Helmholtz free energy is defined as ^> = U — TS, where 
U is the internal energy and S the entropy. The first law of thermodynamics states 
dlJ = 5Q — PdV, where SQ is heat change, the second law for reversible systems 
states 5Q = TdS. Thus = dU - TdS - SdT = SQ - PdV - TdS - SdT, 
so P = — and $ = — J PdV. Since the isothermal compressibility satisfies 
[3 = -£f£ > 0, we get |^ > 0. From (0.2) we get 

(0.3) * unif = - J PdV = -RT\n(V - b) - 

and the convex envelope of (0.3) gives Maxwell's equal area rule for the isotherms 
of (0.2). However, this and related statistical mechanics arguments lack a proof 
of phase transition as defined above. To be more precise, suppose that Maxwell's 
envelope is on the line segment [Vi, V2]. Why does the system take on only V\ and 
V2, but not other values inside the segment? In other words, the implicit meaning 
of arguing this way is that assuming there is a phase transition, it is determined by 
the equal area rule. 

As is often the case in mathematics, we can get a proof of phase transition with 
Maxwell's equal area rule by studying a related larger structure. This was done in 
[20, p. 40]. In a related work, the author also obtained 4/3 as the critical exponent 
[21, p. 360]. Other methods have reported 3/2 (van der Waals), 1.38 [28] and 
(1.21,1.32) [42]. 

We now improve the derivation of van der Waals. 

Let p = h denote the density. The free energy of a nonuniform state p(x) of a 
nonreversible process is a functional of the Helmholtz form I(p) = U(p) — TS(p). 
The cumulative free energy is: 



(0.4) 



^nonunif ^ = f pf u«if^ = f [_ RTp \ n ( p -l _ b ) _ ap 2] 

Ja Jo 



U(p) cannot be only Jq(C — ap)p, but rather — J(x — y)p(x)p(y)dxdy ~ 

Jo P^P — a P — ^-p")dx. Since the term Cp can be neglected due to the mass 
constraint, one can conclude that 



(0.5) 



I{p) = * nonun if (p) - J 1 -^p"p + Jj-RTpHp- 1 -b)- ap 2 
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Note that — ap 2 — RTp(ln(p 1 — b)) is concave (double- well) for p and T satisfying 
RT > (<) 2ap(l — pb) 2 . However, note that the nonlocal and the — Jq ap 2 terms ap- 
peared quite separately, therefore it may not be entirely convincing that one is an ex- 
tension of the other. Also, if — J Q ap 2 is replaced by — J Q J Q J(x — y)p(x)p(y)dxdy, 
the solutions of the EL equation are always continuous, so we have no phase tran- 
sition as defined above. 

In his Nobel lecture, van der Waals expressed his absolute conviction that molecules 
associate in complexes not of chemical origin. He called them pseudoassociations. 
Thus the derivation can be improved by adding the nonlocal term to (0.4): 

(0.6) r nod (p) = - f f J(x - y)p{x)p{y)dxdy + C \—RTp ln(p _1 - b) - ap 2 ]. 
Jo Jo Jo 

The mass constraint makes possible an addition of a linear term, which we choose 
so that — RTplnlp^ 1 — b) — ap 2 has equal depth wells. Also, on a bounded in- 
terval boundary effects may come into play, so J is not necessarily translationally 
invariant: J = J(x,y). Now (0.6) is qualitatively the same as derived in [4] 

(0.7) I(v) = ±J J J(x,y)(u(x)-u(y)) 2 dxdy + J W(u(x))dx, 
where W is the double-well function 

(0.8) W{u) = -^ju 2 - i« 2 + kT[(l + u) ln(l + u) + (1 - u) ln(l - u)], 

with j(x) = f J(x 1 y)dy and k the Boltzmann's constant. For some works on this 
and related models see the References, especially the collective diffusion kinetics 
of the phase transition in polymer gels, what the authors in [24] consider one of 
the most exciting problems in current condensed matter physics, and a nonlocal in 
time evolution discussed in [37]. 

Let G(u) = -\u 2 + kT[(l + u) ln(l + u) + (1 - u) ln(l - it)], G* be the convex 
envelope of G and 

(0.9) I*{u) = -\[ J[u](x)u(x)dx+ f G*(u(x))dx. 

2 Jo Jo 

Let g* — G*', [u,u] be the interval on which g* is constant and v* = g*{u) = 
g*(u). We define 



Let it be a local minimizer of (0.7). Let v = g{u) and xo be such that v(xq) = v*. 
We now show that under some conditions v'(xo) ^ 0, which implies that the set of 
discontinuities of u is finite. 

Let So > be such that J(x,y) > ^J(x,x) for all (x,y) such that \x — y\ < 25o- 
Let S E (0, 5 ) and Is = [x — S,xo + S\. We define 

v + = 5 2 + max{w(a;)|a; e Is}, v~ = mm{v(x)\x E Is} — S 2 . 

Let 

(0.10) <t>={u- s{v-))xi 5 - YZY 8 (J i ( u -s( v ~))x[o,i]\is- 
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Since J Q <p — 0, we have I(u) < I(u — (j)), or 

> J {G*(u) - G*(u -<(>)- 4>J[u] + \<t>J[4>]}- 

On Is we have 

G*(u) -G*{u-(j>) -v<j>> (g*(u -(j>)- v)<f> > (v~ - v+)<f> 

and 

i ^ m] > j(^o) ( ^ ^ _ # 

On [0, 1] \ /a we have 

G*(«) - G*(u - <j>) - v<j> > (g*(u - 0) - s'W)^^ / 

and 

2 W/ ^ J 2(1-2*) + 2(1-2^ 

Taking into account these inequalities we obtain 

j- J(xo,x ) /[o,i]\jj /[o,i]\/ 6 J _ max^gj, j(x) _ I[o,i]\i 5 , 
~ I 4 2(l-2(5) 2 1-25 (l-2<5) 2 J ^ 



Denote the expression in the curly brackets by C(J,g,5). After dividing by fj 
we get 

v+ -v~ > C(J,g,8) [ (u- s(v-)) + 0(5 2 ). 
Jls 



In a similar manner, we take 

= (u-s(v+))xi s - Y^25^jj ( u "^ w+ ))^[o,i]\J. 



Us 

and obtain 

v + - v- > C( J, g, 8) [ (s(v+) -u) + 0{5 2 ). 
Jls 

Adding these two inequalities and dividing by 28 we get 

V+ ~ V ~ > g^% ,+) -s(v-)) + 0(6). 



28 

Letting 8 — > we get 



\v'(x )\>^^-(u-u), 



where 

C(J, g) = + £li _ j(a;o ) _ jf 1 ^ (tt) . 

Thus v'(xo) ^ if C(J,g) > 0. For now this shows that any local minimizcr does 
not take values in (u, u) and proves the existence of phase transition. 



PERIODIC MICRO-STRUCTURE 



5 



The well-balanced scaling 
(O.H) J £ (*,y)-ij s (^) -eJ l (x,y), 

proposed in [31], where J s > and J R \x\ J s {x)dx < oo, can be justified as repre- 
senting an interplay of attractive (J s ) and repulsive (J 1 ) chemical forces. This is 
consistent with an established view in physical chemistry, e.g., "But the establish- 
ment of a well-defined periodicity between the lamellae, or the rods, depends on the 
existence of long-range forces between them: attractive (van der Waals) or repulsive 
(electrostatic, steric, plus the short-range Marcelja repulsion)" [12, p. 2296]. 

Its qualitative shape of that of the "mexican hat" , a notion that comes from 
mathematical biology, where such kernels are also often used. However, it is in- 
teresting to observe that since we assume only J R \x\J s (x)dx < oo, J e can be 
nonnegative. It is then its shape, not sign, that gives pattern formation. 

The construction of periodic minimizers with discontinuous interfaces now fol- 
lows from a sequence of lemmas. Let I f denote / with J f given by (0.11). Then 
e~ 1 I e T-converges as e — > to a singular limit defined on L 2 (0, 1) by 

I (u) = l ^■^^-\titij\x,y)(u(x)~u(y)Ydxdy u £ BV((0, 1), {-1, 1}), 
\ oo otherwise, 

where c = inf {M: „ (±oo)=±1} J R J R J s (x - y)(u(x) - u(y)) 2 dxdy + J R W(u(x))dx and 
1 1 Dm 1 1 is the variation measure of u [1]. Here T-convergence is denned as 

1. For every {u e } C L 2 (0, 1) with lim e ^o u e = u, liminf e ^o £ Ie(u e ) > -^o(w); 

2. For every u £ L 2 (0, 1) n BV((0, 1), {— 1, 1}), there exists a family 

{u e } C L 2 (0, 1) such that lime^oWe = u , and limsup e ^ £ _1 / £ (m £ ) < Io(u). 

These two inequalities often go together the compactness property 

3. Let e n be a sequence of positive numbers converging to 0, and {u n } a 
sequence in L 2 (0, 1). If e~ 1 I en (u n ) is bounded above in n, then {u n } 
is relatively compact in L 2 (0, 1) and its cluster points belong to 
W((0,1),{-1,1}). 

A further related property is that a strict local minimum u of I perturbs to a local 
minimum u e of e _1 / e . This was shown and used in [23] to obtain local minimizers 
on dumbbell domains. In the proof, by a standard argument, a minimizer u e of 
e _1 / e is first constructed in a small closed ball around uo. Then using properties 
1 — 3 it is shown that for small enough e > 0, u e lies in the interior of the ball, thus 
is a local minimizer of e~ 1 I e in L 2 . 

When G is nonconvex this argument goes through a convexification and the 
calculation above. It is shown that a local minimum of Iq in a subspace of 
BV((0, 1), {— 1, 1}) having a fixed number of jumps is also a local minimum of 
Iq in L 2 (0, 1). This observation reduces the problem of determining local minima 
of /g to a finite dimensional one. Since |(||Dti||(0, 1)) is equal to the number of 
jumps of u, it is enough to investigate only the second term of Iq. In general, its 
critical points are determined from a system of nonlinear algebraic equations and 
this is not an easy, or indeed, pleasant task. 

The Ohta-Kawasaki functional follows from [2], though it has been shown in 
[27] that it is not sophisticated enough to predict the complicated phase separation 
morphologies. However, to confirm that pseudoassociations exist and to calculate 
the critical point exponent, we make a study in one dimension. 
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The robust functional 

(0.12) I(w)= [ [\e 2 w"{xf + W{w'{x))+w 2 {x)]dx 1 

Jo 2 

can model e.g., the twinned martensite phase in Nitinol [25]. The elastic foundation 
third term in (0.12) is nonlocal. Namely, let us consider (0.12) with the boundary 
conditions w(0) = w(l) = 0. Let w = v', v = (— L> 2 ) _1 (u — m), where 



-D 2 : {v e W 2 - 2 : v'(0) = w'(l) = 0, f 1 v = 0}^{weL 2 : [ 

Jo Jo 



w = 0}. 



Then w' = v" — m — u, w" = —u' and 

fw 2 = f 1 v 12 = - f 1 v"v = [\u-m){-D 2 )- 1 (u-m)= C m(— £> 2 ) _1 u 
Jo Jo Jo Jo Jo 

and (0.12) becomes 

I(u) = [ [\-€ 2 u'(x) 2 + W(m - u(x))]dx + f f G(x,y)u(x)u(y)dxdy, 
Jo 2 J J 

with G as in [32]. Using the improved van der Waals derivation, we obtain (0.7) 
with J = J e - the strong separation diblock copolymer functional. With G and 
small e - asymptotically periodic local minimizers with sharp interfaces. 



References 



G. Alberti and G. Bellettini, A non-local anisotropic model for phase transitions: asymp- 
totic behaviour of rescaled energies, European J. Appl. Math. 9 (1998), 261-284. 
S. Alexander and J. McTague, Should all crystals be bec? Landau theory of solidication 
and crystal nucleation, Phys. Rev. Lett. 41 (1978), 702. 

J. M. Ball and C. Mora-Corral, A variational model allowing both smooth and sharp 
phase boundaries in solids, Commun. Pure Appl. Anal., in press. 

P. W. Bates and A. Chmaj, An intcgrodiffcrential model for phase transitions: stationary 
solutions in higher space dimensions, J. Statist. Phys. 95 (1999), 1119-1139. 
P. W. Bates, P. C. Fife, X. Ren and X. Wang, Traveling waves in a convolution model 
for phase transitions, Arch. Ration. Mech. Anal. 138 (1997), 105-136. 

F. Bethuel, G. Orlandi and D. Smets, Convergence of the parabolic Ginzburg-Landau 
equation to motion by mean curvature, Ann. of Math. (2) 163 (2006), 37-163. 
M. Bode, Spin-polarized scanning tunneling microscopy, Rep. Prog. Phys. 66 (2003), 523- 
582. 

D. Brandon and R. C. Rogers, The coercivity paradox and nonlocal fcrromagnetism, Con- 
tin. Mech. Thermodyn. 4 (1992), 1-21. 

J. Carr and R. Pego, Invariant manifolds for metastablc patterns in ut = t 2 u X x — f(u), 
Proc. Roy. Soc. Edinburgh Sect. A 116 (1990), 133-160. 

X. Chen and Y. Oshita, An application of the modular function in nonlocal variational 
problems, Arch. Rational Mech. Anal. 186 (2007), 109-132. 

R. Choksi and P. Sternberg, On the first and second variations of a nonlocal isopcrimetric 
problem, J. Reinc Angew. Math. 611 (2007), 75-108. 

P. G. De Gennes and C. Taupin, Microcmulsions and the flexibility of oil/water interfaces, 
J. Phys. Chcm. 86 (1982), 2294-2304. 

J. de la Figuera, F. Leonard, N. C. Bartelt, R. Stumpf and K. F. McCarty, Nanoscale 
periodicity in stripe-forming systems at high temperature: Au/W(110), Phys. Rev. Lett. 100, 
186102 (2008). 

M. del Pino, M. Kowalczyk and J. Wei, On De Giorgi conjecture in dimension N > 9. 

H. F. Ding, W. Wulfhekel and J. Kirschner, Ultra sharp domain walls in the closure 
domain pattern of Co(0001), Europhys. Lett. 57 (2002), 100-106. 

[16] Y.-F. Duan, H.-P. Liu, L. Yl, Atomically sharp domain walls of Rb 2 Cd 2 (504)3 langbeinites: 
Ad electron effects, Physics Letters A 350 (2006), 278-282. 



[7. 

[»: 

[9 

[io: 
in 
[12: 

[13 



[14 
[15 



PERIODIC MICRO-STRUCTURE 



7 



[17] J. L. Ericksen, Equilibrium of bars, J. Elasticity 5 (1975), 191-201. 

[18] R. L. FOSDICK AND D. E. MASON, On a model of nonlocal continuum mechanics. II. Structure, 
asymptotics, and computations, J. Elasticity 48 (1997), 51-100. 

[19] G. FuSCO AND J. K. Hale, Slow-motion manifolds, dormant instability, and singular pertur- 
bations, J. Dynam. Differential Equations 1 (1989), 75-94. 

[20] S. Janeczko, A note on singular submanifolds, J. Geom. Phys. 2 (1985), 33-59. 

[21] S. Janeczko, T. Mostowski and J. Komorowski, Phase transitions in ferromagncts and 
singularities, Rep. Math. Phys., 21 (1985), 357-381. 

[22] G. D. Kahl, Generalization of the Maxwell criterion for van der Waals equation, Phys. Rev. 
155 (1967), 78. 

[23] R. V. Kohn and P. Sternberg, Local minimizers and singular perturbations, Proc. Roy. 

Soc. Edinburgh Sect. A 111 (1989), 69-84. 
[24] E. S. MATSUO and T. Tanaka, Kinetics of discontinuous volume-phase transition of gels, J. 

Chcm. Phys. 89 (1988), 1695. 
[25] S. Muller, Singular perturbations as a selection criterion for periodic minimizing sequences, 

Calc. Var. Partial Differential Equations 1 (1993), 169-204. 
[26] K. Nakamura, Y. Takeda, T. Akiyama and T. Ito, A. J. Freeman, Atomically sharp 

magnetics domain wall in thin film Fe(110): a first principles noncollincar magnetism study, 

Phys. Rev. Lett. 93, 057202 (2004), 
[27] P. D. Olmsted and S. T. Milner, Strong-segregation theory of bicontinuous phases in block 

copolymers, Phys. Rev. Lett. 72 (1994), 936. 
[28] A. PAROLA and L. Reatto, Hierarchical reference theory of fluids and the critical point, 

Phys. Rev. A 31, (1985), 3309. 
[29] G. B. Partridge, Wenhui Li, Y. A. Liao and R. G. Hulet, M. Haque and H. T. C. 

Stoof, Deformation of a trapped fermi gas with unequal spin populations, Phys. Rev. Lett. 

97, 190407 (2006). 

[30] M. Pratzer and H. J. Elmers, M. Bode, O. Pietzsch, A. Kubetzka and R. Wiesen- 
DANGER, Atomic-scale magnetic domain walls in quasi-onc-dimcnsional Fe Nanostripes, Phys. 
Rev. Lett., 127201 87 (2001). 

[31] X. Ren and L. TruSKINOVSKY, Fine scale microstructurcs in 1-D elasticity, Journal of Elas- 
ticity 59 (2000), 319-355. 

[32] X. Ren and J. Wei, Wriggled lamellar solutions and their stability in the diblock copolymer 
problem, SIAM J. Math. Anal. 37, 455-489. 

[33] P. ROSENAU, Free-energy functionals at the high-gradient limit, Phys. Rev. A 41 (1990), 
2227. 

[34] O. Savin, Regularity of flat level sets in phase transitions, Ann. of Math. 169 (2009). 

[35] S. Serfaty and I. TlCE, Lorcntz space estimates for the Ginzburg-Landau energy, J. Func. 
Anal. 254 (2008), 773-825. 

[36] J. Serrin, The area rule for simple fluid phase transitions, J. Elasticity 90 (2008), 129-159. 

[37] L. TruSKINOVSKY and A. Vainchtein, Quasicontinuum models of dynamics phase transi- 
tions, Contin. Mech. Thermodyn. 18 (2006), 1-21. 

[38] M. S. Turner, M. Rubinstein and C. M. Marques, Surface-induced lamellar ordering in a 
hexagonal phase of diblock copolymers, Macromolcculcs 27 (1994), 4986-4992. 

[39] R. van Gastel, N. C. Bartelt, P. J. Feibelman, F. Leonard and G. L. Kellogg, Re- 
lationship between domain-boundary free energy and the temperature dependence of stress- 
domain patterns of Pb on Cu(lll), Phys. Rev. B 70, 245413 (2004). 

[40] G. M. Whitesides and M. Boncheva, Beyond molecules: Self-assembly of mesoscopic and 
macroscopic components, Proc. Natl. Acad. Sci, USA 99 (2002), 4769-4774. 

[41] W. Windl, T. Liang, S. Lopatin and G. Duscher, Modeling and characterization of atom- 
ically sharp "perfect" Ge/Si02 interfaces. Materials Science and Engineering B 114-115 
(2004), 156-161. 

[42] R. A. YOUNG, Theory of quantum-mechanical effects on the thermodynamic properties of 
Lcnnard-Joncs fluids, Phys. Rev. A 23 (1981), 1498. 

Division of Integral Equations, Faculty of Mathematics and Information Science, 
Warsaw University of Technology, Pl. Politechniki 1, 00-661 Warsaw, Poland 
E-mail address: A . Chmaj (Smini . pw . edu . pl 



